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Recently, much work has been devoted to the calculation of O(a) corrections to the decay rate 
of pionium, the -k + tt~ bound state. In previous calculations, nonrelativistic QED corrections were 
neglected since they start at order a 2 in hydrogen and positronium. In this note, we point out that 
there is one correction which is actually of order a times a function of the ratio ii r ct/m e , where fi r 
is the reduced mass of the system. When \i T ct -C m e , this function can be Taylor expanded and 
leads to higher order corrections. When fi r a ~ m e , as is the case in pionium, the function is of 
order one and the correction is of order a. We use an effective field theory approach to calculate 
this correction and find it equal to 0.4298 aTo- We also calculate the corresponding contribution to 
the dimuonium n~ bound state) decay rate and find a result in agreement with calculations by 
Jentschura et al. 



I. INTRODUCTION 

The decay rate of pionium (tt + tt~ bound state) has generated much interest recently, because of the possibility of 
extracting from a future experiment at CERN Jl]] |^] precise values of the pion scattering lengths which can, in turn, 
be related to the quark condensate parameter, a fundamental ingredient of ChPT. In order to do so, however, all 
corrections of order a times the lowest order decay rate must be evaluated. These fall in several categories and have 
been the subject of several recent papers || [§ H ■ 

One can distinguish between pure QED corrections and corrections which arise from the interplay of strong and 
QED effects. The former are associated with an expansion in a whereas the latter contains an expansion in powers 
of a and other ChPT coefficients, such as scattering lengths, the pion mass difference, etc. 

The pure QED corrections can be subdivided in relativistic and nonrelativistic contributions. This separation of 
scales is used, in one form or another, in all nonrelativistic bound state calculations but is the most transparent in the 
effective field theory (eft) formalism. The eft appropriate to the study of nonrelativistic QED bound states, NRQED, 
was developed by Caswell and Lepage more than ten years ago j|] . NRQED has since been used to calculate high order 
corrections to the hyperfine splitting of both muonium and positronium g. Using NRQED, it is easy to convince 
oneself that in hydrogen and positronium, the 0(a) QED corrections are purely relativistic in nature whereas the 
nonrelativistic corrections start at order a 2 . This argument was used in [|j to dismiss the nonrelativistic corrections^. 

In this note, we want to point out that the vacuum polarization correction to the Coulomb interaction (referred to 
as the "VPC" correction in the following) is actually of order a times a function of the dimensionless ratio fj, r a/m e , 
where fi r is the reduced mass of the bound state. In hydrogen and positronium, this ratio is very small (~ a and a/2, 
respectively) and the function can be Taylor expanded, leading to a correction of higher order. In pionium, however, 
the ratio is equal to 0.997 and the function turns out to be of order one. This means that the VPC correction 
represents an 0(a) correction to the lowest order rate. We want to emphasize that the observation that the VPC 
correction is enhanced in atoms with large reduced mass is not new; it is, for example, discussed in the Landau and 
Lifshitz textbook on quantum electrodynamics in the context of energy levels. 

Another system in which this type of correction is large is dimuonium, the /i + fi~ bound state. In that system, 
the ratio n r a/m e is about 0.754 and the VPC correction is also of order a times the lowest order decay rate. This 
contribution was evaluated, among several other corrections, in a recent paper by Jentschura et al Jxcf] . In this note, 
we use effective field theories to calculate the VPC contributions to both dimuonium and pionium. Our dimuonium 
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result, 0.337aro, agrees with pQ] Our result for pionium is equal to 0.430 oTq, which is as large as the other corrections 
evaluated in § [§. 

We first consider the VPC corrections to the decay rate of dimuonium. For the present calculation, the only relevant 
interactions of the NRQED Lagrangian are 



"-f ermions = ^ + ^ e ^o)j"! 



J 2 

£f ermions = ^ T { idt + — eA ) }tp + same terms with ip 

c 4 ^aa 2 X* ■ X i cr 2 o i ^- c 5 ^ a 2 X* ' X 1 ^ (1) 



where ip refers to the muon field, x refers to the antiunion field and m n is the muon mass. Our sign convention for 
the coefficients c 4 and C5 differs from [Tlfl in order to avoid minus signs in the Feynman rules. To the order of interest, 
we can approximate the gauge derivative D by p since the correction due to transverse photons are suppressed (see 
|jl2|| ). The Feynman rules of the relevant interactions are presented in Fig.[l]. We work in Coulomb gauge in which 
the Coulomb photon propagator is simply given by —l/k 2 . The vacuum correction to the Coulomb propagator is also 
illustrated in Fig.[l]. We consider only the VPC because it is the only interaction which contributes at order a times 
a function of fx r a/m e . The vacuum polarization correction to the transverse propagator can be neglected because 
the vertices connecting transverse photons to fermions are suppressed by at least one power of the muon mass, which 
translates into additional factors of a in the bound state calculation (the interested reader is referred to [0 for more 
details on NRQED power counting rules). 

The four-fermion operator proportional to C4 leads to the decay rate of ortho dimuonium (dimuonium with no 
angular momentum and a total spin equal to 1) whereas the other four- Fermi operator contributes to the decay of 
paradimuonium (L = 0, S = 1). In a nonrelativistic system, the Coulomb interaction is nonperturbative and must be 
summed up to all orders, which is equivalent to solving the Schrodinger equation. We will focus on the decay rate of 
the ground state (n = 1), for which the wavefunction is given by 



87r l/2 7 5/2 

*® = -^P)®^®V ( 2 ) 



where the £'s are the two component Pauli spinors of the muon and antimuon and 7 = fi r a represents the typical 
bound state momentum (the ground state energy is given by — 7 2 /(2// r )). 

In NRQED, bound state properties are computed in two steps. First, one fixes the NQRED coefficients by matching 
QED and NRQED scattering amplitudes, and then one computes the bound state energy by applying time ordered 
perturbation theory to the NRQED vertices (with the Coulomb interaction defining the unperturbed problem). The 
decay rate is given in terms of the imaginary part of the energy by T = -2 Im (E) . 

To illustrate the use of NRQED, we now focus on the decay rate of orthodimuonium. To simplify the following 
discussion, we define the operator 

0( 3 5i) = Im(c 4 ),S 2 (3) 

where the spectroscopic notation 2S+1 Lj is used to indicate the quantum numbers of a state annihilated by this 
operator. In terms of this operator, the decay rate is given by 

Tfi Ct 

r(l 3 5x) = -2 (*|G( 3 5 1 )|*) = -4|*(r = 0)| 2 Im( C4 ) = - Im(ci) (4) 

where a factor of 2 comes from the spin average of S 2 and \^(f = 0)| 2 = 7 3 /7r. The imaginary part of the coefficient c 4 
contains an infinite expansion in powers of a, the lowest order contribution being found by performing the matching 
illustrated in Fig. [2]. All diagrams are computed exactly at threshold (no external momenta), in the center of 
mass frame of the particles, and a nonrelativistic normalization is used for the QED diagram (i.e. a factor of 
1/V2E = 1/V2m is provided for each external particle). Also, we match the amplitude matrix M. and not — iM.. 
The superscript 1 on c 4 indicates that this is a contribution coming from a one-loop QED diagram. The NRQED 
diagram is simply 2Im(c 4 1 ^) whereas the imaginary part of the QED diagram is -2vra 2 /(3m 2 ). Solving for the 
NRQED coefficient, we get 

Mci 1 ^-^ (5) 



3m 2 



which, after putting in Eq.(||) gives the lowest order rate, 
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(0) 



We now turn to the contribution from the vacuum polarization correction to the Coulomb photon. The first thing to 
do is to revise the matching of Im(c4) beyond one-loop. At first, this might seem unnecessary since we are interested 
in 0(aTo) corrections only, but we must be careful about factors of m M /m e which might compensate for factors of 
a. Indeed, this situation occurs in the matching illustrated in Fig. [3]. The result for the imaginary part of the QED 
diagram can be found from ||l3f : 

2 x ( " )2 WW + 4 m, + 22 rn.\ W 
7r \ 16 m e 3 m e 36 m M / 3 

where the term in curly brackets is simply the lowest order result and the overall factor of 2 takes into account the two 
permutations of the diagram. The term proportional to m p /m e represents a correction equal to 0.283 a and should 
therefore be treated as an 0(a) correction instead of an a 2 correction. To the precision of interest in this work, we can 
neglect the log term, which comes from the running of the coupling constant, and all the other terms not enhanced 
by a factor of m M /m e . The calculation of the NRQED one-loop scattering diagram in Fig. [3] is straightforward: 



2?r a 2 , r d 3 k -2^ r f -4aV(l - v 2 /3) 3 m p 2 , 2tt a 2 
'~3~m^' J (2tt) 3 p J V k 2 {l-v 2 )+Am 2 e _ 16^T a ^~~3~^f 



which is seen to cancel exactly the m M /m e term in the QED diagram. This is not a coincidence. The m^/ffie term is 
associated to the low energy behavior of the QED diagram, and since NRQED is designed to reproduce QED in the 
nonrelativistic limit, this term had to be present in the NRQED diagram. This will hold for higher order loop diagrams 
as well. The end result is that the coefficient Im(c4) does not receive any correction of order a(am fM /m e ) n ~ a from 
the matching to QED. Therefore, we can simply use the lowest order coefficient hn(c^) given in Eq.(^|) in the bound 
state calculation. 

The bound state contribution of the VPC interaction is found by applying second order perturbation theory: 

ST vpc ( 3 S 1 ) = 2(n\0( 3 S 1 ) *V I") (9) 

where V vpc is the potential corresponding to the vacuum polarization correction to the Coulomb interaction: 

f 1 , -4aV(l-w 2 /3) 
Jo fc 2 (l - t> 2 ) + 4to 2 

For the sum over intermediate states, we must include the full Coulomb Green's function G(p, q) in terms of which 
the previous expression can be written as 



flVrSi) = -8x 



d 3 p' 8^/ 2 7 5 / 2 -Tra 2 
(2tt) 3 (p* 2 + 7 2 ) 3m 2 
d 3 pd 3 qd 3 k 

(frr)» ° M 

1 , -Aa 2 v 2 (l-v 2 /'S)\ 8^/2 7 5/2 

dv— = (11) 

£; 2 (1 - v 2 ) + Am 2 e ' (P~ k ) 2 + 7 2 



= 2 / ^y* ^ s) ( t dv -y v2{i - v2/3) ) »r r (12) 

J (2tt) 9 \J k 2 {l - v 2 ) + Am 2 J (p- k) 2 + 7 2 

where, in the first expression, a factor of -2 comes from the relation T = —2 Im(E), a factor of 2 comes from the spin 
average, and another one comes from the two sides on which the interaction can take place. The second expression 
gives the correction in terms of the lowest order rate, Tq. For the Coulomb Green's function, we use the following 
expression derived by Schwinger pif : 

G(p, q) = J^(2n) 3 6 3 (q) - — 4 R(p, q) (13) 
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where the first term corresponds to no Coulomb interaction in the intermediate state, the second term corresponds 
to one Coulomb line and the third term corresponds to the the exchange of two or more Coulomb photons. We will 
refer to these three terms as the — C, 1 — C and R contributions, respectively. A graphical representation of Eq. ( |l2] ) 
is presented in FigJ4l. An explicit expression for R(p, q) can be found in [T^ ]. All we need for the present work is the 
following integral |lq] 

[(d 3 P d 3 q)f(p)R(p,q) = -4ttV r d P X^L (H2) ~ I + ^Arctan(g) - iln(l + ^) + (14) 
J Jo (p +7 ) V 2 P 7 2 Y P +TJ 

where /(p) is an arbitary function of the magnitude of the three- momentum p. 

We now calculate the correction to the decay rate of orthodimuonium. The no Coulomb correction to Eq.(|l2|) can 
be integrated analytically and is found to be 

AEo-c = 4 Jo(r) «r (15) 

7T 

where Iq is a function of the dimensionless ratio r = "//m e = [x r a/m e which is 0.754 in dimuonium. For the sake of 
completeness, we give Io(r) for arbitrary values of r: 



,/ 5 1 tt,1 1 * (-2r + l/r + l/r 3 ) f A 
h = ir { n + - (t + ^-o) + , — ^ <^ ArcTan 



TV 



^(-32+15tt) (r = l) 



(r<l) 



= 7T 2 



5 1 7T/1 Is (-2r + l/r + l/r 3 ) f , V / ^ _ T + r N l\ , 



Using Eq.(^5|) with r = 0.754, the no Coulomb contribution to the dimuonium decay rate is found to be 

A^o-c (dimuonium) = 0.216ar . (17) 

We have reduced the one Coulomb and R term contributions of Eq.([l2]) to two dimensional integrals (using the 
identity (|l4|)) which were evaluated using the adaptive Monte-Carlo integration routine Vegas jl6). The results are 
respectively 

A£ , i_ c (dimuonium) = 0.08674 aT (18) 



AE R (dimuonium) = 0.0344 ar (19) 

leading to a total equal to (sum of Eqs.([l5|), ( |l8[ ) and (|l9|)) 0.337aTo. This result agrees with Ref. [^| which gives 
(1.06/7r)aro = 0.337oTo. Notice that the lowest order decay rate of paradimuonium is also given in terms of a four- 
Fermi interaction (the interaction proportional to C5 in the Lagrangian) so that the VPC correction is also 0.337 oTq 
where Tq is now the lowest order paradimuonium decay rate. 

The calculation of the VPC correction to the pionium decay rate is very similar to the dimuonium calculation. This 
is the most apparent by using a nonrelativistic field theory reproducing ChPT at low energy, which we will refer to 
as "NRChPT" . The only NRChPT interactions relevant to the present calculation are the Coulomb interaction (for 
the charged pions) and the four-pion vertex of the form 7r + 7r~7r 7r (a more detailed presentation of NRChPT and 
calculations of the other 0(a) corrections will be presented elsewhere). 

The matching of the four pions interaction is presented in Fig. [5]. At threshold, the ChPT amplitude for the 
7r + 7r~7r°7r° vertex is equal to 327t(ciq — Oq)/^' wnere a o is the scattering amplitude in the isospin channel i. Using a 
nonrelativistic normalization for the external states, we divide by a factor 

(y/2E^) 2 y/2E„ +y /2E„- w Am^m^ w 4m^ + (20) 

where we have set Am^ = m w + — m w o =0. As we will see below, this is possible because we are interested in the 
decay reate to leading order in Am v (the 0(Am 7r ) corrections will be treated in a future publication). The end result 
of the matching is that the coefficient of the NRChPT four-pion interaction is 

8^ (qg - qg) 

^4- pions — 9 I -^-L ) 

6 m + 
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For the coulomb photon and the VPC interaction, we use the same Feynman rules as in the dimuonium calculation. 
As before, the Coulomb interaction must be summed up to all orders, which means that the tt + tt~ state will be 
described by the same ground state Schrodinger wavefunction as before. The lowest order rate is trivially computed 
by considering the NRChPT interaction sandwiched between wavef unctions: 



AE = \j Wf m Cipions l Wf G ^° m Cipions l Wf W) (22) 

where the overall factor of 1/2 is a symmetry factor and G^o^o (q 2 ) is the free propagator of the 7r°7r° pair, given by 

G A »(?*) = w 5——- (23) 

Eq - E intermedlat e 2to^+ - 2m,o - q 2 /m v o + le 

where we have neglected the binding energy in Eq (the full expression is Eq — 2m 7r +(l — a 2 / 8)) since it leads to a 
correction of order a 2 To. The difference 2vn^ — 2ra^a = 2Am 7! appearing in Eq.(^3|) cannot be set to zero since it 
leads to the imaginary part of the integral and, therefore, to the first order decay rate. A trivial integration gives, for 
the integral Eq.(|22|) 



, Wm^o (a[] — a 2 ,) 2 



AE = - ^(O)! 2 ^" V "° to4 "° ; (A + t|V*W2A^) (24) 

where A is an ultraviolet cutoff on the \q\ integration which would be canceled by a counterterm in a calculation of the 
bound state energy fll7[ . Since we are only interested in the decay rate, we need only to concentrate on the imaginary 
part of the energy. Working in first order in Am,, Eq.(p3) finally leads to a decay rate equal to the well-known result 



To = -2IME) = 1^/ ^(08 - 08)' a 
9 y m w + 

Turning now to the VPC correction, we have to carry out the integral given in Eq. (p"2|) , using the reduced mass of 
pionium, /i r = /2. The no-Coulomb term is clearly given by Eq.([l5|), with r — m T +a/(2m e ) ~ 0.997. We find a 
result equal to 

A£ -c (pionium) = 0.26677 ar . (26) 

We have also computed the one Coulomb and R term contributions to the pionium decay rate using VEGAS and 
found 

AEx^c (pionium) = 0.1193 aT (27) 

AE R (pionium) = 0.04372 qT . (28) 
Our final result is the sum of Eqs.@, (|27]) and @: 

AE vpc (pionium) = 0.4298 aT (29) 
which is of the same order of magnitude as the other corrections calculated in fl and 0] . 
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FIG. 1. NRQED Feynman rules (q represents the charge of the particle). The last three rules can also be applied to NRChPT. 
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FIG. 3. Vacuum polarization correction to the coefficient Im(c4) (the imaginary part of the coefficients is implicit in the 
figures) . 
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FIG. 4. VPC correction to the dimuonium decay rate. The three diagrams represent the no Coulomb, one Coulomb and R 
contributions (the imaginary part of the coefficients is implicit in the figures). 
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FIG. 5. Lowest order matching of the ChPT and NRChPT 7r + 7r~7r°7r° interaction. 
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